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A polychromatic microwave quadrature interferometer has been characterized using several laboratory plasmas. 
Reflections between the transmitter and the receiver have been observed, and the effects of including reflection 
terms in the data reduction equation have been examined. An error analysis which includes the reflections, 
modulation of the scene beam amplitude by the plasma, and simultaneous measurements at two frequencies has been 
applied to the empirical database, and the results are summarized. For reflection amplitudes around 10%, the 
reflection terms were found to reduce the calculated error bars for electron density measurements by about a factor 
of 2. The impact of amplitude modulation is also quantified. In the complete analysis, the mean error bar for high- 
density measurements is 7.5%, and the mean phase shift error for low-density measurements is 1.2°. 
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I. Introduction. 

A millimeter wave interferometer which uses a 
beam consisting of three discrete frequencies has been 
used to measure phase shifts in laboratory plasmas. This 
instrument was designed and built as a primary 
diagnostic tool for plasma experiments in support of 
research on propulsion technologies. The design was 
derived from several sources, notably including the 
UCLA Large Plasma Device 70 GHz quadrature 
instrument [1], a 35 GHz (heterodyne) interferometer 
used at Forschungszentrum fUr Mikrostrukturtechnik 
[2], and discussion of early microwave interferometer 
experiments by Heald and Wharton [3]. The microwave 
frequencies provide good sensitivity, and incorporation 
of phase shifters and additional detectors provides the 
benefits of quadrature capability [4]. Generally, 
interferometry at multiple frequencies is useful in 
species discrimination [5], multiple fringe identification, 
and error analysis. In the present instrument, the three 
frequency components, while spatially overlapped, are 
essentially independent, and the interferometer may be 
used for mono- or bi-chromatic measurements without 
operational modification. Also, although not addressed 
here, the interferometer has reflectometry capability, 
which, along with the polychromaticity, is suited to 
operation near critical density, as in, for example, 
investigation of collision frequencies in transient 
plasmas [3]. 

The scene beam of this instrument is intended 
for applications in which the transmitter and receiver are 
separated by a distance of approximately one meter. 
Depending in part on the alignment, it was observed that 
reflections set up between the transmitter and receiver 
resulted in significant effects on the phase shift 
measurements [3]. Probable candidate causes for this 
circumstance are thought to be (1) the polychromatic 
design of the interferometer, in which all three 


frequencies travel in a common wave guide, (2) the 
optical configuration, which is essentially symmetric 
with respect to transmitter and receiver, and (3) the 
flexibility in alignment due to the use of cables, rather 
than waveguides, in connecting the transmitter and 
receiver. The inclusion of these reflections in the data 
reduction model was found to reduce the errors 
calculated for the phase shift measurements. 

An analysis of data collected with this 
interferometer is presented here, based largely on the 
results of electron density measurements made in two 
experimental devices, a high-density magnetized plasma 
column and a low-density electrostatic discharge. The 
high-density plasma is the core of the Variable Specific 
Impulse Magnetoplasma Rocket (VASIMR) propulsion 
testbed and is generated in a long-pulse mode (a few 
sec.) [6]. The low-density device is a spherical Inertial 
Electrostatic Confinement (IEC) chamber operating in a 
sustained discharge mode [7]. A total of 3685 phase 
shift measurements are included in the database, taken 
from 312 experimental runs in plasmas from three 
species. The phase shifts measured ranged from < 0.5° 
to 530°, corresponding to areal electron densities from 
= 10 u cm' 2 to =5xl0 13 cm" 2 . These phase shifts 
were determined with the standard interferometer 
equation as well as with an equation including 
reflections, and the results are compared. For the present 
analysis the VASIMR data are evaluated at two 
simultaneous frequencies, and the IEC data at two 
frequencies sequentially. Also, attenuation of the scene 
beam, or more generally, modulation of its amplitude, 
during plasma operations is important in many cases. A 
Monte Carlo error analysis was performed for this data, 
and statistical distribution functions for the error bars, as 
well as for the phase shift differences among models, 
are compiled from the results. 
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II. The Instrument 


The interferometer consists of a transmitter and 
receiver connected by coaxial cables. The microwaves 
are generated on the transmitter and reference signals 
for the three frequencies are transmitted, at reduced 
frequencies, over the connecting cables to the receiver, 
where they are multiplied before combination with the 
scene components. The working frequencies are 70, 90, 
and 110 GHz, and the scene beam spatially superposes 
these frequencies in a nominally Gaussian transverse 
profile using Millitech Gaussian focusing and collection 
optics. The separation of the transmitter and receiver is 
variable about the nominal 1 meter and the scene beam 
waist is located approximately 60 cm from the 
transmitter. The full-width at half-maximum of the 
waist is about 1.8 cm at 70 GHz. For each of the three 
frequencies, the reference and scene components are 
each split into two signals before combination at the 
receiver, and one of the two is shifted in phase to 
produce a quadrature arrangement. A schematic 
drawing is shown in Fig. 1. 

The spatial distribution of microwave irradiance in 
the scene beam was mapped for the 70 GHz component 
using an Epsilon Lambda crystal detector. The detector 
face is a 2.9 cm diameter metal disk with a 0.27 x 0.55 
cm rectangular aperture in its center. The detector was 
placed in front of the transmitter optic and translated 
axially and transversely while recording voltages output 



Figure 1 Schematic of the microwave interferometer. 
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Figure 2 Measured spatial distribution of microwave 
irradiance for the 70 GHz beam. The measurement grid 
spacing is 1 cm axially and 0.16 cm transversely. 

from a preamp attached to the detector. The transverse 
translation was accomplished using a stepper motor and 
the axial translation was performed manually. In the 
mapping, a standing wave pattern was observed 
between the transmitter and the detector, giving rise to 
axial oscillations of the measured irradiance of the 
microwaves at the detector. Fig. 2 gives the results of 
the mapping. The horizontal axes give the axial and 
transverse positions, in cm, and the vertical axis gives 
the signal strength (in mV). The axial position is 
measured from the transmitter optic, and the transverse 
position from an arbitrary zero point chosen near the 
transverse maximum. The axial spatial resolution of 
Fig. 2 is 1 cm and the oscillation due to the reflection is 
not resolved, but is evident at an aliased period of s 2.9 
cm (A/2 = 0.214 cm). Fig. 3a gives a plot of a single 
transverse profile taken at an axial distance of 60.2 cm, 
and Fig. 3b gives a plot of the axial profile near the 
waist, scanned with a stepper motor at the central 
transverse position. The axial oscillations of Fig. 3b 
have an amplitude approximately 30% the mean signal. 

During normal operation of the interferometer, 
reflections similar to those observed in the beam 
mapping, but less marked, were also observed. These 
reflections occur between the transmitter and receiver, 
and their strength, and the degree to which they affected 
the interferometry, varied with the alignment of the 
instrument and with the experimental configuration. 
Also, the reflections were usually less apparent in the 
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FIG. 3. Measured transverse (a) and axial (b) profiles of the 
70 GHz beam. The discrete points are the measurements. 
The data in (a) are from Fig. 2, at an axial position of 60.2 
cm. The oscillations in (b) are due to standing waves set up 
by the detector. 

90 GHz beam than at the other frequencies. In 
discussing the reflections, it is important to note the 
distinction between spatial oscillation and amplitude 
modulation. The oscillations of Fig. 3b are due, like 
ideal interferometer signals, to the phase relation 
between waves of constant amplitude. In contrast to this 
are changes in the amplitude of the collected scene 
beam signal at the interferometer detector due to, for 
example, changes in the source output or modified 
coupling due to the plasma [3], 

The interferometer is equipped with a rotary 
phase shifter which varies the optical path length, but 
not the physical path length, for the scene beam 
upstream of the transmitter optics. However, because 
the plasma itself (or other scene target) varies the phase 
of the scene beam downstream of the transmitter optics, 
calibration of the interferometer was performed by 
translating the transmitter along the optical axis z with 
a stepper motor, called here the phase motor [1]. In 


particular, the effect of scene beam reflections, which 
depends on the optical path length between the 
transmitter and receiver, is in this way incorporated in 
the calibration. The translation was made possible 
because of the flexible connection between transmitter 
and receiver. The general procedure is to perform a 
calibration scan just prior to the plasma run. Then the 
transmitter is stationed at the measurement position z s , 
and data are acquired during the formation and 
dissolution of the plasma. The calibration scans were 
usually made over a range from z s -d to z s +d t 

where d is chosen to be slightly larger than half of the 
maximum wavelength. The rotary phase shifter, while 
not used for calibration, was used during alignment and 
it was generally observed that the effect of reflections 
was not evident in the phase shift pattern produced by 
this device. 

One problem we have encountered in the use of 
this instrument is failures of the detector diodes. In the 
course of several tens of hours of operation at three 
different laboratory facilities, damage to these diodes 
resulted in the loss of a frequency channel on four 
occasions. The manufacturer (MilliTech) has suggested 
that these failures are related to high-voltage pulses 
which occur episodically and originate in plasma 
support hardware, such as power supplies. Replacement 
of the diodes is comparatively straightforward, but in 
practice it entails a change of the quadrature shift for the 
affected channel. In addition, there is a degree of 
variability in the quadrature shift, particularly for the 70 
GHz channel, which is observed in the absence of diode 
failures. For both of these reasons, the receiver was 
eventually modified by the addition of adjustable 
shifters, so that the quadrature may easily be set in 
general. However, this adjustment was not available for 
the data presented here, and phase shifts as small as 30° 
and as large as 116° have been used in consequence. 

II. Experimental Plasmas 

This interferometer has been used to make 
measurements in several experimental plasmas. A 
database taken from two of those experiments is used 
for the work here. The first of these is the VASIMR 
VX-10 device (at JSC in Houston TX), which creates a 
plasma column by flowing gas down a quartz 
confinement tube through a helicon antenna. The 
plasma column is shaped by solenoidal magnetic fields 
and has a length of approximately 2 meters. A 
schematic is given in Fig. 4. The data used here were 
part of a series of experiments on diffusion and 
recombination in the plasma conducted over a period of 
, four days, and were collected at an axial position at the 
exit of the confinement tube [6]. The ion-cyclotron 
resonance heating (ICRH) antenna, normally installed 
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just downstream of this position, was removed for these 
experiments (it is not shown in the figure). The plasma 
column diameter at this position was approximately 2.4 
cm, and the gas was He with typical flow rates around 
230 seem. Most of the data were taken at the centerline 
of the column, but several transverse scans are also 
included. The VASIMR device was operated in a pulsed 
mode and the typical duration of a “shot” was 
approximately 2.3 seconds. 

The second experimental device was a spherical 
IEC chamber, which creates a plasma by static 
discharge in a neutral gas and confines the plasma using 
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Figure 4 Schematic of VASIMR microwave interferometry. 
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Figure 5 The low density IEC plasma device showing the 
interferometer calibration and plasma scanning 
arrangement (top). 


a radial electrostatic field [7]. A sketch of the IEC 
device is given in Fig. (5). The vacuum chamber is 60 
cm in diameter, and is jacketed for temperature 
regulation, producing plasmas which are stable over 
extended periods of time (typical experiments last many 
10’s of minutes). It is equipped with two 180° viewports 
having clear apertures of 10 cm diameter. The data used 
here were collected using two different central cathodes, 
15 and 23 cm in diameter, respectively, over a range of 
transverse locations within the viewport access. Typical 
values for the confinement voltages were - 1 kV, and 
for the chamber pressures = 1 Pa. Two gas species, Ar 
and D 2 , are included. 

III. Data Analysis 

In the phase shift analysis here, it is assumed that 
each of the three frequencies is isolated from the other 
two. The single frequency which is treated explicitly 
may thus represent any of the three. For this single 
frequency, the input channel will be indicated by 
default, and primes will be used to denote the 
quadrature channel. Underlining denotes plasma 
quantities. 

The electric field at the microwave detector is a 
superposition: 

E = E r e~ im +aE s e ii *+* a ~ m) [\+a 2 p e ' (2 ^ P) J (1) 
where E and co are the field amplitude and frequency, 
respectively; the subscripts R and S denote the 
reference and scene fields, respectively; a and p are 
the amplitude modulation and reflection coefficients, 
respectively; <p is the phase shift of the scene beam 

relative to its phase at the position z = 0; fa is the 
phase difference between the reference beam and the 
non-reflected scene beam at the zero position; and fi is 
the phase shift of the reflected component of the scene 
beam relative to the non-reflected component at the zero 
position. The amplitude coefficient a is included to 
account for cases in which there is modulation of the 
scene beam amplitude when the plasma is present (e.g., 
refractive bending); by definition a = 1 when the 

plasma is absent. Letting E —> y jyce Q E , where y is the 
detector responsivity (in volts watt" 1 m 2 ), c is the speed 
of light, and e 0 is the permittivity of free space, the 
voltage output from an interferometer detector is thus 
modeled as: 

V = E 2 +(.aE s ) 2 ^ + {a 2 p) 2 ]+V D + 

2aE R E s [cos(0+^ o ) + a 2 /9 cos (30 + (^ + /?)] (2) 

+ 2a 4 pE 2 cos(20 + P) 

where V D is the “dark” voltage. For a quadrature 
channel, Eq.(2) is unchanged except that the zero phase 
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0 O flo + 0q = 0Q * where is the quadrature shift. 
For p = 0, Eq. (2) reduces to the standard form [1]. 

Eq. (2) is the basis for the phase shift 
measurements. The constants, etc ** needed to 

evaluate the phase in Eq. (2) are denoted Qu , and are 
treated, for convenience, as the components of a vector 
C 0 , with basis c k , k = l,K. Other related quantities, 
identified with the subscript k , will also be treated here 
as vector components in this basis. The C 0Jc are listed 

in Table 1 for the three cases used in this work; /r is a 
scale factor included in the phase term for the 
calibration, </> = kz . The separate reflections case, in 
which the input and quadrature channels are treated 
independently, is not physically distinct from the 
coupled case, but is computationally faster. The 

theoretical voltage V , then, is a function of C 0 and the 
phase <j>\ for the calibration data, V is treated as a 
function of the phase motor position z , 

V = V(z;C 0 ) (3a) 

and for the plasma data, as a function of <p = KZs + A# 

V = V(</>;C Q ) (3b) 

where A<t> is the plasma phase shift 

The calibration constants C 0 ,C' 0 were determined 
from the calibration scans by least squares fitting of the 
scan data to the model. Specifically, the amplitude 
residues R, R’ for the fitting are defined as 



n=l 


where V n ,V' are the microwave signal voltages 
measured at the phase motor positions z„ for n = l,N , 

and the C 0 ,C' 0 are chosen to minimize /?,/?'. Examples 
of calibration scans and the curve fits based on Eq. (4) 
are given in Fig. 6 for each of the three frequencies. 
These examples were chosen specifically to illustrate 
the instrument reflections. The reflection coefficient p 


is defined as the ratio of the electric field amplitude, 
incident at the detector, of the reflected component to 
that of the unreflected component, and the average 
value for the cases shown in Fig. 6 is p - 8.6%. The 

fractional RMS residue Rq=(R/ N) U2 /(AV 12) , 
where AV is the peak-to-peak voltage of the sinusoidal 
oscillation of the calibration scan, has average values of 
6.2% for the standard model and 1.3% for the reflection 
model. Again, the deformations caused by p > 0 were 
not observed in phase scans performed with the 
upstream rotary shifter. 

For the plasma measurements, the scene phase is 
determined by minimization of the local residues 

r+r'= [v s -V(*;C 0 ) ] 2 + fo' -V(frC' 0 ) f . (5) 

In this case, the constants C 0 ,Cq are fixed to their 
calibration values and the minimization is accomplished 
either by (1) the variation of 0 alone, or (2) by 
simultaneous variation of 0 and a . The plasma phase 


shift is then given by the difference, Atp = 0 - 0 , (6) 
for paired voltages with the plasma ( V S ,V^ S ) and 
without the plasma (V S ,V^). For the case p - 0, an 

analytic solution ^(v) exists [1]. 

Monte Carlo simulation was used to estimate the 
uncertainty in the phase shift measurements. The 
contributing error sources were taken to be uncertainties 
in (1) the calibration constants, and (2) the plasma 
measurement voltages. The calibration constant 
uncertainties represent systematic errors, “fossilized” 
from the calibration data [9]. Each of the fitted 
calibration constants C 0tk was taken to be the mean of 


Gaussian distribution f(C k ) with standard deviation 
<T k . A theoretical residue R was defined. 


tf(O s £[y(vC 0 )-v(z,,;C)] 2 

=i 

dV l 2 

w t iC '- c A 



(7) 


and the a k were determined from the conditions that 
(1) the mean value of R be equal to the measured 


Table 1 Calibration constants for fitting of calibration scans to the model, Eq. (2), for 3 cases. 


Reflections 

No. of 
parameters 

Qu 


None 

2x4 

Er + Es 

E' R 2 + E? + V^2E' R E' s ,ti,y 

Separate 

2x7 

E'R > E$ > 05 > K y V D , /?, P 

E'^E's^yy^p',? 

Coupled 

11 

E R yE Sy <j\ i ,Ky D ,p, P 

E' R ,E' s ,<&-^Ky D ,p,p 
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Figure 6 Calibration scans for the interferometer. The circles 
are the input channel voltages V and the diamonds are the 
quadrature channel V f . The broken lines are fits to the 
standard model and the solid lines are fits to the reflection 
model. The reflection coefficients p , and fractional RMS fit 

residues Rq>R^ for the standard and reflection cases 
respectively, are: 

70 GHz: p = 12.6%, R 0 =8.5, 1.7%, =9.7, 1.9%; 

90 GHz: p = // = 5.6%, Rq = 3.7,0.71%, /?£ = 4.2, 0.66%; 
110 GHz: p = 7.7%, R 0 =6.0, 1.7%, R^ =5.1, 1.4%. 


residue, 


K 

l 

k~\ n~\ 


N ^v ^ 2 


dc k 


= R 


( 8 ) 


and (2), the individual errors be proportional, 
statistically, to ( dV /3 C k )~‘ , 



For the measurement voltage, the standard deviation 
u v of the distribution f(V) about V s was estimated 
using both the scatter of steady state data and the local 
residues: 

<7v=V y iL + r, (10) 

where v^. is an RMS value of V s (r), taken directly 
from measurement data for the IEC case and estimated 
from pre/post test data for the VASIMR case (see Sec. 
V). Distributions were then calculated from 

Q.i s Cojc ± V2 y t a k ^ 

V t mV s ±j2y,0 v 

where the >’, are derived from random numbers X i 

between 0 and 1 such that erf(y ( ) e ~‘ 2 dt = x i , 

and the index i is for the statistical sampling. From 
these distributions, distributions for the phases <p and 0 

were calculated and fit to Gaussian curves. The error bar 
for the phase shift measurement is then given by the 
root-sum-square (RSS) of the fitted widths, 

SA0 = . (12) 


The electron density is related to the phase shift by 



where k x = 2n! A is the wave number, n(z) is the 


density and n A is the critical density, at which the 
plasma becomes opaque [8]. For n(z) « n A , the phase 
shift is proportional to the line integral of the density 
and does not depend on the path length or the shape of 
the density distribution along the path. Eq. (13) derives 
from the refractive index of the plasma and assumes 
v c « (o , where v c is the effective collision frequency, 
which is expected to be the case for the data here (e.g., 
VASIMR p ~ 0.1 Pa, T e ~ 5 eV). Also, the microwave 
polarization for the VASIMR measurements was 
parallel to the axial magnetic field, so that, in the cold 
plasma model (thermal velocity v, h «v^~c ), the 

dispersion is not affected by this field (= 0.4 T at the 
measurement position, giving co\la> 2 = 0.026 at 70 
GHz, where co B is the gyrofrequency) [3]. 

Eq. (13), along with the phase shift errors, Eq. 
(12), then, gives electron densities and density errors for 
each of the frequency channels of the interferometer. 
For the analysis here, two simultaneous frequencies are 
used for the density measurements. For a single data 
point, the two frequency channels give error 
distributions /p/ 2 , with means Mi’l-h an< ^ widths 
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S ( ,S 2 » which may be combined as a scalar product 
/,/ 2 to obtain a reduced distribution with mean 



Formally, this may be motivated by interpretation of the 
densities in the individual distributions as paired 
measurements, one from each frequency channel. If the 
two distributions have the same mean and width, for 
example, the reduced distribution has also the same 

mean and a width reduced by -J2 ; that is, the reduced 
width scales as a standard deviation of means for N = 2 
[9]. For the general case, in which Mi* M 2 an( * 
Si * S 2 , a local bias error is defined such that 

P\ + P2 = |m -/' 12 I+K " Mn\ . subject to 0 < # <* S,- , 
and the widths of the individual distributions become 
of =S?+$ , i = 1,2 , (15) 

giving the reduced deviation cr 0 , 

% 2 = a i 2 + °2 2 ’ ( 16a ) 

and mean 

Mo=Mn- (16b) 

Thus this correction fi t accounts for any discrepancy in 
the means Mu M 2 while preserving Mn * obtained 
directly from the original distributions. The reduced 
distribution is then integral-normalized, 

jfi(Mi>^ux)f 2 (M 2 ^ 2 ’ x ) dx = e ~ r ' 2 > (17a) 

where 



vanishes for Mi = Mi- The integral, Eq. (17), which 
essentially gives the overlap of the two distributions, is 
used below in comparing the final results among the 
models. 

In the following section, these Gaussian 

distributions / , associated with individual 
measurements, will be termed populations, to 

distinguish them from the distributions of population 
widths to be discussed. 

IV. Results and Discussion 

The VASIMR database represents essentially all 
the helium shots performed during the test series at the 
confinement tube exit position. It is comprised of 102 
separate shots, including a total of 3475 individual 
phase shift measurements. The measurement sample 
rate was approximately 13 Hz for most of the data, but 
higher rates, up to = 23 Hz, were used in some cases. 


For the analysis here, the intent was to provide a 
representative statistical sample, and only coarse 
screening of the data for anomalies was performed. 
Ideally, a calibration scan would be performed for every 
plasma shot. In the current system, however, constraints 
on the communications speed for recording the position 
of the transmitter along the optical axis limited the rate 
of calibrations, so that the ratio of calibration files to 
data files for this data is = 0.6. These constraints also 
limited the resolution for the calibration scans; a typical 
resolution was 0.012 cm, corresponding to 10° at 70 
GHz. 

Post-test, the calibration scans were fitted to the 
model, Eq. (2), using both the standard case and the 
coupled reflections case (Table 1) in Eq. (4). Each shot 
file was then reduced using the calibration constants 
from both models at each data point in the file, and this 
reduction was performed twice, once with fixed a - l 
and once with (X determined from the best fit to Eq. (5). 
A Monte Carlo simulation using 5000 samples was then 
performed for each of the 4 results at a point, giving a 
set of 4 error bars at each point. The error bar is defined 
here as the unit width (half-width at e~ l maximum = 

•J2a ) of the Gaussian population of phase shifts. 

The reduced phase shift data were then used to 
compile distribution functions. All the data from all the 
shots were pooled together and parsed into intervals, or 
bins, based on the value of the phase shift or on the 
value of the amplitude modulation coefficient. The 
modulated-amplitude reflection model is taken as the 
reference case, and values from this model are used for 
the intervals, for both A<p and a . Then the mean value 
of the error bar in each interval was computed. Also, the 
mean phase shift differences between each of the 
models and the reference case were computed for each 
interval. Generally, the change in amplitude, as well as 
the magnitude of reflections, is larger for the low 
frequency channel. For reference, mean and extreme 
values for the reflection and amplitude coefficients in 
this database are given in Table 2 (the value a = 0.05 is 
an arbitrary minimum set in the data reduction, and 
occurred identically only a few times). The mean 
quadrature shifts for the two frequencies were 
approximately 31 and 116°, respectively, for the 70 and 
110 GHz data. 


Table 2 Mean and extreme values for the reflection and 
amplitude coefficients for the VASIMR data. 



p 

a 

70 GHz 

110 GHz 

70 GHz 

110 GHz 

Min 

0.103 

0.0577 

0.05 

0.54 ] 

Max 

0.191 

0.110 

1.27 

1.28 

Mean 

0.131 

0.0782 

0.498 

0.934 
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The distributions of the error bars over phase shift 
are shown in Fig. 7. The means for the 70 GHz data are 
14.2 and 14.7°, respectively, for fixed and modulated 
amplitude in the standard model, and 9.3 and 6.8°, 
respectively, in the reflection model. For the 1 10 GHz 
channel, the corresponding numbers are 8.1, 6.5, 5.3 and 
3.0°, respectively. The average number of 
measurements per bin is about 54. Fig. 8 shows the 
distributions of phase shift differences. The differences 
show some phase dependence which is common to both 
frequencies; for example, they are comparatively small 
around -130°, and comparatively large near -200°. The 
large oscillation evident in the 70 GHz data is due in 
part to the attenuation, which is more substantial at the 
low frequency (see Table 2). Essentially all the 70 GHz 
data with \A<p\ > 320° have a< 0.5, and omission of 
this data also substantially reduces the difference 
extrema for | A</> | < 320°. The large differences for the 
70 GHz data are also due in part to the departure from 
quadrature. The error bar distributions are fairly uniform 



Figure 7 Distribution of the phase shift error bars for the 
VASIMR data. The squares are fixed amplitude, the circles 
modulated amplitude. Open symbols are the standard model, 
solid symbols include reflection. The error bar values are unit 
widths (half-width at 1/e maximum). 


between frequencies, and show a mild phase 
dependence qualitatively similar to that of the 
differences. However, high frequency maxima in the 
error distributions are notably absent in the variable 
amplitude models. The distributions over amplitude 
coefficient, not given in the figures, showed no 
pronounced correlations; there was a weak tendency for 
the errors to decrease with increasing a< 1, particularly 
for the standard models. 

For comparison. Table 3 gives the results of the 
data reduction and error analysis applied to several 
small sets of synthetic data. A fixed set of phases <p ^ , 
giving shifts from 0 to 281°, was used to generate 
perturbed phase data <p e , <p' e by adding to 0^ random 
phase errors from a Gaussian population with a unit 
width of = 5°. From the <p £y <j>' e synthetic 

voltages V y V' were calculated for 8 shot files of 84 
points each, using four different quadrature shifts for the 
two cases 1 and 0.66. Calibration files were 




Figure 8 Distribution of the phase shift differences, relative to 
the reference model, for the VASIMR data. The squares are 
fixed amplitude, the circles modulated amplitude. Open 
symbols are the standard model, solid Symbols include 
reflection. Lines with symbols are mean differences and lines 
without symbols are RMS differences. 


8 






Table 3 Mean values, in degrees, of phase discrepancy magnitudes \<t> - | , in italics , and phase shift error bars SA<p , in 

bold, for reduction of synthetic data with ~ 0.1. The true values are \<f> € -0^) = n~ xn 5 $ ^ = 2.82° and SA<p^ - 
2 W2 ^tnic = ™7°. 


Quadrature 
angle </> Q 

^true 

Standard (p- 0) 

Reflections (p> 0) 

Fixed a 

Modulated a 

Fixed a 

Modulated a 

30° 

1.0 

3.92 

113 

5.08 

11.8 

1.98 

835 

3.25 

8.28 

0.66 

20.7 

12.0 

3.96 

13.6 

16.6 

10.1 

3.43 

11.2 

57° 

1.0 

3.93 

104 

4.65 

103 

2.15 

7.81 

2.56 

7.69 

0.66 

15.4 

20.6 

3.22 

113 

13.6 

21.1 

2.58 

9.78 

90° 

1.0 

4.24 

10.9 

4.60 

10.4 

2.21 

8.11 

2.28 

8.07 

0.66 

8.17 

23.0 

3.37 

11.0 

8.32 

20.8 

2.51 

9.19 

120° 

1.0 

4.39 

123 

4.68 

12.5 

2.24 

8.54 

2.32 

8.45 

0.66 

14.6 

19.0 

3.25 

12.7 

14.2 

15.6 

2.38 

9.96 


synthesized in a similar manner using the same %>*• 
The resulting data were then processed to produce a set 
of “measured” phases <p , phase shifts A<j > , and error 
bars SA<p for each of the 8 shot files. The mean phase 
error, as determined by , has values of 9.41° 

(fixed- a) and 4.10° (mod- a) for the standard models 
and 7.66° (fixed- a) and 2.66° (mod- a) for the 
reflection models, which compares with the true value 
of 2.82°. The mean phase shift error bar SA<j> for the 
reference case is 9.08°, as compared with the true value 
7.07°. Substantial variation in the results due to the 
departure from quadrature is generally confined to 
fixed- a reduction in the attenuation cases. In these 

cases, shows a phase dependence which is 

similar to that of the difference distributions in Fig. 8; 
the magnitude of the effect increases with departure 
from quadrature and is at a minimum for <j> Q = 90°. 

In calculating electron densities from the VASIMR 
data, the objective here was to determine the fractional 
error bars and differences. The details of the spatial 
distribution of the plasma along the optical axis will 

impact the absolute densities for n of order n x , but 

have little effect on the fractional errors. For the present 
purpose, then, the plasma distribution is modeled as a 
step function with a width of 2.45 cm. This diameter is 
the equivalent width for a Gaussian distribution 
determined from inversion (uncorrected for the beam 
profile) of transverse scans made with the 
interferometer, and is consistent with the visual 
appearance of the plasma. Using this step function in 
Eqs. (12) and (13) gives for each phase shift 
measurement a density and a density error for each of 
the two frequencies, 70 and 110 GHz. These densities 
were normalized to the low-frequency critical density. 


n >l( 70 GH 2 ) =6.1xl0 13 cm 3 , and used to compile 

distribution functions, just as with the phase shifts. The 
results for the error bars for combination of the two 
frequencies as per Eq. (16) are shown in Fig. 9. As a 
check on the error analysis. Fig. 9 also gives the overlap 
integral Eq. (17) for the each of the three models 

compared with the reference model, that is, e~ yn 
with /, and f 2 as reduced density error populations for 
the comparison and reference models, respectively. The 
mean values of e~ r ' 2 are 0.80 for the fixed amplitude 
models and 0.95 for the modulated amplitude. In 
general, the density error bars are reduced by inclusion 
of reflection and modulation. Table IV gives mean 
values for the data of Fig. 9, along with mean values for 
the fractional density differences for the two 
frequencies, where the fractional difference is defined as 

An = (/I, iogh* - «70GHz ^” 1 10GHz + n 70GHz )] an( * * S 

independent of the error analysis. The behavior of 
| An | for the different models is qualitatively similar 
to that of the error bar Sn , although the effect of 
modulation is more pronounced than that of reflection 
for the differences. An itself indicates the tendency 

n 110GHz <n 70GHz- 

For the EEC data, a single phase shift measurement 
was based on data collected during several contiguous 
plasma cycles, each having a typical duration of 
approximately 10 seconds. Such measurements were 
made over a range of neutral background pressures and 
transverse locations, and separate calibration files were 
obtained for each measurement. The voltages from the 
interferometer (V S ,V^ S and V S ,V $ ) used for phase 

determinations were averages taken over steady state 
windows within the plasma cycles. Given the voltages 
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FIG 9 Distributions for the combined frequency (70 /1 10 
GHz) VASIMR electron density measurements, (a) The 
normalized error bar Sn / k^oghz) * w ^ ere <5 n is the unit 

width, (b) The integral of the product of error populations, 
Eq. (17), for comparison with the reference model. The 
squares are fixed amplitude, the circles modulated 
amplitude. Open symbols are the standard model, solid 
symbols include reflection. 


Table 4 IV Mean values, in percent, for the fractional error 
bars and differences for the VASIMR density data at 70 
and 110 GHz , from the data of Fig. 9(a). 



Standard (p= 0) 

Reflections(/?> 0) 

Fixed 

a 

Modulated 

a 

Fixed 

a 

Modulated 

a 

combined 
frequency 
error bar 

19.6 

14.4 

10.2 

7.54 

difference 

between 

frequencies 

-0.076 

-2.11 

-3.49 

-4.92 

magnitude 

of 

difference 

22.9 

16.2 

17.9 

14.3 


and the calibration data, then, the data reduction 
proceeds just as with the previous case. The total 
database is comprised of measurements made over a 


period of several months, and only measurements made 
at transverse positions > 5 cm from the viewport center 
were excluded. At 5 cm, which is the nominal edge of 
the window, the scene beam signals are reduced by 
clipping typically to 10 - 15% of their central values, as 
determined from the calibration scans. The total number 
of phase shift measurements in the EEC database is 210. 
The measurements are monochromatic but data at two 
frequencies, 70 and 90 GHz, were used in separate test 
series. The quadrature shift was 54° for the 90 GHz data 
and varied from 24 to 39° for the 70 GHz data. The 
calibration resolution was typically about 0.0034 cm (3° 
at 70 GHz), higher than for the high-density data, and 
these scans were fit to the separate reflections model 
(Table I). 

Histograms of this data are presented in Fig. 10. 
Most of the shifts are less than 2°, and all are less than 
6°. In this case the difference distributions were 
compiled using the average of the results for the four 
models as the reference phase. The distribution means 




Figure 10 Distributions of phase shift measurements 
A(j> (heavy line, left axis), phase shift differences zl0-d0 ave 
(broken lines, right axis), and error bars SA(j> (solid lines, right 
axis), for the IEC data, taken in Ar and D2. 
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Table V. Mean differences and errors, in degrees, for the IEC 
i phase shift data at 70 and 90 GHz. 




Standard (p = 0) 

ReflectionsGo > 0) 

Fixed 

a 

Modulated 

a 

Fixed 

a 

Modulated 

a 

70 

GHz 

magnitude 

of 

difference 

0.599 

0.556 

0.564 

0.553 

error bar 

2.65 

2.57 

1.42 

1.35 

90 

GHz 

magnitude 

of 

difference 

0.0988 

0.0931 

0.0890 

0.0902 

error bar 

1.71 

1.52 

1.26 

1.03 


are given in Table V. It is seen that the difference 
distributions show little variation, and that the error bars 
are fairly insensitive to amplitude modulation, which 
was very small in this data. The mean values of the 
reflection coefficient were 0.032 and 0.019, 
respectively, for 70 and 90 GHz. The fractional change 
in the calculated error bar with reflection terms remains 
significant, but the absolute size of the error is small in 
this data, less than 3° for all cases. 

VI. Conclusion 

The cable coupling of the transmitter and 
receiver in this instrument permits flexibility in the 
experimental configuration and, particularly, simple 
translational phase calibration. The translational 
calibration, in turn, permits modeling of residual 
reflections which occur between the transmitter and 
receiver and results in reduction of the measurement 
uncertainty. The quadrature configuration, in addition to 
eliminating sensitivity nulls and resolving phase change 
direction ambiguities at extrema, is critical in 
determining (monochromatic) phase shifts in the 
presence of significant amplitude modulation, which is 
typical in high-density data. The polychromaticity of the 
interferometer has proved useful here in error analysis 
for electron density measurements, by permitting 
combination of error distributions for simultaneous 
measurements at two frequencies, and in providing 
independent support for the analysis, from comparison 
of the simultaneous measurements themselves. 

In the VASIMR experiment, temporally 
resolved interferometer data from plasma shots were 
analyzed to determine phase shifts and error bars using 
the standard phase model and a modified version 
including reflection terms. For both models, the data 
were also processed neglecting amplitude modulation 
by the plasma, for comparison. The measured phase 
shifts are from 0 to 530°, and the differences in phase 
shift among the models are typically 5 to 20°, with 
larger differences occurring for strong attenuation off 
quadrature. The phase shift error bars, calculated from a 


Monte Carlo algorithm, show a weak phase dependence 
and have means from 3 to 15°; in the full model the 
mean errors are 6.8° and 3.0°, respectively, for the two 
frequencies, 70 and 110 GHz. The inclusion of 
reflections reduced the average phase shift errors by a 
factor of 1.8 and the inclusion of modulation by a factor 
of 1.3. The error bars for the low frequency are larger 
than those for the high frequency by approximately a 
factor of 2, which compares to the (asymptotic) ratio 1.6 
for the phase shifts themselves. Electron densities were 
calculated from the phase shifts at the two frequencies, 
and the difference in density between the two 
frequencies has mean magnitudes of 23% for the 
standard phase model and 14% for the full model. 
Combining the two frequency measurements in the 
uncertainty analysis results in mean density error bars of 
20% and 7.5%, for these respective models, with a 
general dependence on the model similar to that of the 
phase shifts. 

In addition to the high-density work, a similar 
phase analysis of a collection of 70 and 90 GHz data 
from the steady-state EEC device was performed. These 
data show very small phase shifts (<2° typical) and 
negligible attenuation, with less variation in results 
among the models. The mean error bars for the full 
reduction were 1.4 and 1.0° for the 70 and 90 GHz, 
respectively. 
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